Sound propagation in a cylindrical Bose-condensed gas 
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We study the normal modes of a cylindrical Bose condensate at T — using the linearized time- 
dependent Gross-Pitaevskii equation in the Thomas- Fermi limit. These modes are relevant to the 
recent observation of pulse propagation in long, cigar-shaped traps. We fin d that pulses generated 
in a cylindrical condensate propagate with little spread at a speed c = yj gn/m, where n is the 
average density of the condensate over its cross-sectional area. 



PACS numbers: 03.75.Fi, 67.40. Hf, 67.57.Jj 

In a recent paper, Andrews et al. ]l| presented results 
for the propagation of sound through the Bose conden- 
sate in an axially-symmctric, cigar-shaped trap. These 
observations were made possible by the large aspect ra- 
tio of the anisotropic trap and the localized generation 
of pulses by a laser beam focused at the mid-point along 
the length of the trap. By imaging the propagation of the 
pulse, they were able to determine the sound speed as a 
function of the condensate density. The measured sound 
speed is approximately equal to the Bogoliubov phonon 
velocity as determined by the maximum density in the 
cloud. The authors Jlj] sketch a qualitative theoretical 
argument in support of this observation. 

The purpose of the present paper is to provide a more 
detailed analysis of the dependence of the sound speed 
on the inhomogeneous distribution of atoms in the trap. 
In the propagation of a pulse for short times, the ends 
of the trapped condensate do not come into play and 
we can therefore simulate the experimental situation by 
considering an idealized cylindrical trap which is uniform 
in the z-direction and which has a harmonic confining 
potential of the form 



V(p) = -mu) 2 p 2 



(1) 



in the radial direction. Here, u>o is the trap frequency. 
Treating the condensate in the Thomas-Fermi (TF) ap- 
proximation , the equilibrium density is given by 



no(p) 



2.9 



(2) 



where g = iirah 2 jra is the interaction parameter and R 
is the radius of the cylindrical condensate related to the 
chemical potential, p = V(p) + gno(p), by p = ^muj^R 2 . 

The dynamics of the condensate will be based on the 
linearized time-dependent Gross-Pitaevskii equation J^] 
in the TF limit B . This equation can be recast as a pair 
of quantum hydrodynamic equations 



-V • (n v) 



dSn 

~dt 
dv 

m— - = —gVSn — W5U , 
at 



(3) 
(4) 



where SU (r, t) represents an externally imposed poten- 
tial. Eliminating the velocity from these equations, we 
obtain an equation for the density fluctuation Sn(r, t) 
first derived (for SU = 0) by Stringari H, 



d 2 Sn q 

— — = ^-V • n V<5n 
ot m 



1 



V • (noVSU) 



(5) 



To study the normal modes of the condensate, we set 
SU to zero and look for solutions which correspond to a 
propagating wave of the form 



5n(r,t) = Sn(p)e l(kz - Ult) 



(0) 



We assume that the density fluctuation in the transverse 
directions depends only on the radial variable p and not 
on the azimuthal angle <f). It is clear that more general so- 
lutions having the angular dependence e" n ^ are possible, 
but these will not be considered in this paper. 

We now substitute this form of the solution into (||) 
and, making use of the equilibrium density in (||), obtain 
the equation (SU = 0) 



oj Sn — -uj q • 



k 2 (R 2 



p 2 )5n--{R 2 
P 



o 2s 9Sn 
3P) ^P- 



- (R 2 p 2 



dp 2 



(7) 



where Sn now represents the spatially-dependent ampli- 
tude defined in (g). The allowed solutions of this equa- 
tion define the dispersion relations, to(k), of the various 
modes of the condensate in the cylindrical trap. 

It is now convenient to introduce the dimensionless pa- 
rameters 

uj = — , k = kR 

LOq 

and to define the new independent variable 



2-^-1 

R 2 ' 



-1< x < 1 



(8) 



When expressed in terms of this variable, the density 
fluctuation, now denoted by y(x), is found to satisfy 



1 



d 

dx 



(l-x 2 ) 



<iy_ 

dx 



1-2 Ira 

-ui y 
2 y 



fc 2 (l-x)y = 0, (9) 



which is in the standard Sturm-Liouville form. For fc — 0, 
this equation in fact reduces to the Legendre differential 
equation 



d_ 

dx 



2\ d V 



ax 



+ Xy = 0. 



(10) 



whose solutions are the Legendre polynomials Pi (x) with 
eigenvalues X = 1(1 + 1), I = 0, 1, . . . Thus the modes of 
the condensate can be labeled by the index I and have a 
limiting frequency at k = given by 



Jf(k = 0) = 21(1 + l)u 2 . 



(11) 



The first two frequencies in the sequence, and 2wo, are 
the two m = mode frequencies found by Stringari |4| for 
a cigar-shaped trap in the limit that the axial frequency, 
to z , tends to zero. More generally, the whole sequence in 
( |TT| ) was recently obtained analytically in this same limit 
by Fliesser et al. ||. 

These various solutions have distinct radial eigenfunc- 
tions. For example, for I = 0, y(x) = a constant, and 
the associated density fluctuation has a radially inde- 
pendent amplitude. This mode corresponds to a local 
adiabatic expansion of the condensate brought about by 
a local change in the chemical potential, Sp,. Such a 
change gives rise to a variation in the equilibrium density 
tiq(p\ R) which is spatially independent. In the limit of 
long wavelengths, it costs no energy to move atoms from 
one region of the trap to another and the mode there- 
fore has zero frequency. The other I ^ modes have a 
finite frequency since they involve radial motion of the 
condensate. (Note that according to (Q), the direction 
of the local velocity is normal to constant density sur- 
faces.) For example, for 1=1, the density fluctuation is 
of the form Sn(p) oc (1 — 2 j o 2 /i? 2 ) which has a node at 

p = R/y/2. Unlike the I = mode, all of these modes 
have a number conserving density fluctuation in the sense 

that f dp pSn(p) = 0. The I = 1 mode at k = is a 
radial breathing mode. 

For k ^ 0, we expand the density fluctuation as 



i 



(12) 



where the normalization of the Legendre functions is cho- 
sen such that 



dxPi(x)Pi l (x)=5 w 



(13) 



Substituting expansion ( |12[ ) into (|9|), we obtain the sys- 
tem of linear equations 



1 



1 



- 1(1 + 1) - ^k 2 j ai + ifc 2 Mivav = , (14) 



where the symmetric matrix Miy is given by 



M w = / dxPi(x)xP v (x). 



(15) 



Since the matrix Mw does not have diagonal matrix ele- 
ments, it is clear that its effect on the mode eigenvalues 
first occurs in second order, so that to lowest order in k 2 , 

ujf(k) = 21(1 + 1)uj 2 + ^(kRoj f + 0(k i ) . (16) 

In particular, the / = mode is phonon-like with a ve- 
locity 



1 R / gno(0) 



(17) 



where no(0) is the (maximum) condensate density on the 
axis of the trap. If we compare this velocity with the Bo- 
goliubov velocity yj gn/m in a homogeneous gas of den- 
sity n, we see that the effective density in the cylindrical 
trap is half the maximum value, which is just the den- 
sity, n, obtained by averaging (fa) over the condensate 
cross-section. 

It is interesting to note that the matrix problem in jl^ ) 
is equivalent to the problem of a quantum rigid rotor in 
a gravitational field, having the Hamiltonian 



L 2 

H= — + MgR(l- cost 



(18) 



Here, L is the angular momentum operator, / = AIR 2 is 
the moment of inertia of the rotor and the gravitational 
potential energy of the rotor mass M is referenced with 
respect to its lowest position (the polar angle 9 is defined 
with respect to this position). The angular cigenfunc- 
tions of this Hamiltonian can be expanded in spherical 
harmonics Yi m (9, <j)) and in the m — subspace, the prob- 
lem is then identical to (|lj) . Thus the parameter k plays 
the role of a 'gravitational field' which perturbs the mode 
energies from their 'free rotor' values. 

The solution of (jlj) is easily obtained for arbitrary k 
by numerical diagonalization. In Fig. 1 we show some 
of the lowest mode frequencies as a function of wavevec- 
tor. The lowest mode is a sound-like mode and exhibits a 
negative dispersion. It can be seen that the group veloc- 
ity deviates appreciably from its long wavelength limit 
once the wavelength is comparable to the diameter of 
the condensate. The mode-coupling induced by the k 2 
perturbation in ( |l4| ) of course becomes more significant 
with increasing k and has the effect of lowering the sound 
speed. This coupling is associated with the interplay of 
the wave-like modulation of the density along the axis 
and the strong inhomogeneity of the equilibrium density 
in the radial direction. 

We now return to the experimental situation Q which 
involves the propagation of sound pulses rather than con- 
tinuous waves. These pulses are generated by switching 



2 



on a laser beam which repels atoms from its point of ap- 
plication at the center of the cylindrical trap. We shall 
assume this perturbation to be weak and consider the 
linear response of the condensate. The perturbation is 
taken to be 



6U(r,t) = U e- z2/a2 6(t) 



(19) 



that is, a gaussian potential with no radial dependence, 
which is switched on at t = 0. The equation defining the 
density fluctuation is given by (||) and can be solved by 
introducing a Fourier representation of the density and 
external potential: 



/°° dk 
— e lkz Sn(p,k,t) 



SU(z,t) 



dk 
2^ 



e ikz SU(k,t). 



(20) 



with 



5U(k, t) = y^aU e- a2k2 ^e{t) = SU(k)0(t) . (21) 

Taking the Fourier transform of (^]) and using the variable 
x introduced earlier, we obtain the equation 



d 2 y f (u^kRf 

where we have defined the differential operator 

1 



L = -2Wn 



dx dx 



~-k 2 {l-x) 



(22) 



(23) 



This is just the differential operator defining the mode 
eigenfunctions in (ph. These eigcnfunctions satisfy 



(24) 



L(j>m{x) = Un(k)(j>m(x) 
and form an orthonormal set, 

dx(f> m (x)(f> n (x) = S mn . (25) 

They can therefore be used to expand y(x, k, t) as 

y(x,k,t) = J2 b ™{k,t)(t> m (x). (26) 

m 

Substituting this expansion into ( p2| ) and taking the inner 
product of the resulting equation with <f) n (x), we obtain 
the equation for a driven harmonic oscillator 



b n +uJ 2 n b n = f n (t) 



(27) 



where the forcing term on the right hand side is defined 
as 



4.9 

fn{k)6{t). 



(28) 



The solution of ( |27| ) for t > with the boundary condi- 
tions b n (0) = and b n (0) = is 



b n (k,t) = ^ n 2 \ k 2 [1 - cosuj n (k)t] . 



(29) 



We have here explicitly displayed the dependence of the 
various quantities on the wavevector k in the Fourier ex- 
pansion. This essentially completes the solution for the 
density fluctuation Sn(p, z,t). 

To analyze the time evolution of the density pulse, it is 
convenient to consider the average of 5n(p,z,t) over the 
cross-sectional area, A, of the condensate. We therefore 
define 



1 



6n(z,t) = — I dASn(p,z,t) 



/OO JJ, i p 1 

^.e^Y^b^k^-J^dx^ix). (30) 

According to (|T2|), the eigenfunctions <fi n (x) are them- 
selves given by the expansion 



4> n (x) = j24 n) (k)Pi(x), 
i 



(31) 



where the coefficients a|™^ (k) defining the n-th eigenvec- 
tor of dT4| ) have the limiting value a\ n \k = 0) = <5;„. (In 
other words, the eigenfunction <p n {x) reduces to P n {x) in 
the k — > limit.) The orthonormality of the Legendre 
functions thus yields 



5n(z, t) 



Similarly, we find 



dk 
— < 
2tt 



ikz 



fn(k) 



(uJokR) 2 fj n ) 
2V2g 



V2 



\k) 



^4 n) (^n(fc,t). 



(32) 



V3 



a[ n) (k) ) SU(k). (33) 



Thus, for small k, only the n — and n = 1 modes are 
excited appreciably. 

Substituting ( p9| ) together with ([53) into (|32j), the av- 
eraged density fluctuation can be eva uated numerically. 
As an example of the calculation, we present in Fig. 2 
results for the propagation of a pulse. The width of the 
gaussian perturbation in these calculations was chosen 
to be a = 1.5i?, which corresponds approximately to the 
experimental situation. Once the pulse is launched it 
is seen to propagate at the speed c given by (|l~7|). The 
shape of the pulse will change in time as a result of the 
/c-dispersion of the normal modes, but the effect is weak 
since the pulse is mainly made up of the long-wavelength 
n = sound modes. The particles carried off by the pulse 
leave behind a static depression at the origin which cor- 
responds to the new equilibrium shape of the condensate 
in the presence of the applied localized perturbation. 
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These theoretical results for the propagation of a pulse 
are in reasonable agreement with the experimental ob- 
servations . One difference is the stronger wave packet 
dispersion observed experimentally which may be due to 
the actual variation of the condensate density along the 
trap axis. However, a more serious discrepancy concerns 
the speed of propagation. In the experiment it was con- 
cluded that the sound speed was determined by the max- 
imum condensate density according to d ~ yj 'gno(0)/m, 
as opposed to the value we find in (|l7|). However it is 
clear from the published figure that the measured sound 
speed actually lies below d at the higher trap densities. 
At these densities, the theoretical value c would lie be- 
low the measured values by about the same amount that 
d lies above. The difference between our theoretical re- 
sult and experiment could be explained by errors in the 
determination of the trap density on which the theoret- 
ical value is based, but the claimed experimental preci- 
sion (W. Ketterle, private communication) would seem 
to eliminate this possibility. 

On the theoretical side, one can question the validity 
of the TF approximation. For the given trap parame- 
ters E[ and a typical density of n (0) ~ 3 x 10 20 m~ 3 , 
the radius of the cond ensate R, in units of the oscillator 
length I — y^h/mujQ, is approximately 5. This implies 
that the TF approximation is providing a good estimate 
of the ground state density over most of the occupied 
volume. The validity of the TF approximation is further 
supported by the fact that the frequencies of the low- 
lying collective excitations in these traps || are close to 
the TF values Q . This observation is consistent with cal- 
culations based on the Bogoliubov approximation 
which show a reasonably rapid convergence to the TF 
limit with increasing TV, the number of trapped atoms. 
However, the situation with regard to pulse propagation 
is somewhat different in that shorter wavelength excita- 
tions are being probed. Within the hydrodynamic formu- 
lation 0| of the Bogoliubov approximation, an additional 
term appears on the right hand side of (0) which repre- 
sents the fluctuation in the 'quantum pressure' associated 
with the kinetic energy. This term becomes increasingly 
important as the wavelength of the density fluctuation is 
reduced. Unfortunately, it cannot be estimated using the 
TF density since the failure of the TF approximation at 
the edge of the condensate (|Jl^] leads to a divergent re- 
sult. Thus, explicit numerical solutions of the Bogoliubov 
equations are required in order to ascertain the accuracy 
of the TF approximation for the conditions under which 
the pulse propagation experiments were performed. On 
the other hand, at long wavelengths, we expect the TF 
prediction for the sound speed, as given by (|l7|), to be 
valid. 
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FIGURE CAPTIONS 

Fig.l: The dispersion of the lowest modes in a cylindri- 
cal condensate as a function of the wavevector k. 
ujq is the radial trap frequency and R is the TF 
condensate radius. 

Fig. 2: The propagation of a pulse generated by a gaussian 
perturbation applied at the origin at t — 0. The 
lowest curve is at a time t — 4/wo, and the interval 
between each successive curve is At = 4/ujq. 
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